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Abstract 

Our goal is to examine the role of anomalies in the hydrodynamic regime of field theories. We 
employ methods based on gauge/gravity duality to examine R-charge anomalies in the hydrody- 
namic regime of strongly t'Hooft coupled, large N , J\f = 4: Super Yang-Mills. We use a single 
particle spectrum treatment based on the familiar "level-crossing" picture of chiral anomalies to 
investigate thermalized massless QED. In each case we work in the presence of a homogeneous 
background magnetic field, and find the same result. Regardless of whether a particular current 
is anomalously non-conserved or not, as long as it participates in an anomalous 3-pt. correlator, 
its constitutive relation receives a new term: j'^ oc —dP'^'^'Bi^ p"^ . This agrees with results found by 
Alekseev et.al. for QED. We include a general, symmetry based argument for the presence of such 
terms, and use linear response theory to determine their coefficients in a model with anomalous 
global charges. This last method, we apply to briefly examine baryon transport in chiral QCD in 
a strong magnetic field. 
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I. INTRODUCTION 

Hydrodynamics jll] provides a universal description of the long range, long time-scale 
behavior of a wide variety of thermal systems. The hydrodynamic quantities are those for 
which small perturbations from equilibrium have relaxation times which diverge as their 
wavelength diverges. As a result, hydrodynamics has been formulated as a classical effec- 
tive field theory with various fields, typically the conserved quantities, playing the role of 
fundamental degrees of freedom. Hydrodynamic behavior has been seen to arise from finite 
temperature quantum field theory 

At the same time however, one feature that sets aside relativistic quantum field theories is 
the existence of quantum anomalies, which cause some classically conserved quantities to be 
non-conserved. An unsolved problem is how to incorporate the effects of quantum anoma- 
lies into the hydrodynamics description of a thermal field theory which contain anomalies. 
In a QCD plasma, the hydrodynamic theory should reproduce the three-point correlation 
functions including the anomalous part (an example of such correlation functions is that 
of the axial vector, vector and baryon currents). Another example is the massless QED 
plasma (or the QED plasma at such high temperatures so that the mass of the electron can 
be neglected). In this case, magnetohydrodynamics has to be enlarged to incorporate the 
axial current j^^ = ip'j'^'j^tp. Although this current is not conserved, d^j^^ ~ E ■ B, at large 
distances the total axial charge should change slowly because the conducting plasma cannot 
support a large long-distance electric field E. 

In this paper, we consider a simpler problem where anomalies enter the hydrodynamic 
equations already at the linearized level. Namely, we consider a theory with a set of global 
conserved charges j"- for which the triple correlators contain anomalous contributions. We 
then turn on an external magnetic field coupled to one of the global charges. An example 
of such a system is high-temperature QCD in a very large magnetic field. Now anomalies 
appear already at the level of two-point correlators, and hence should be manifested in the 
linearized hydrodynamic equations. 

We shall argue that in the presence of a background magnetic field coupled to some 
conserved charge, the constitutive equation for the currents is modified in the presence of 
quantum anomalies. In an anomalous theory, any current j'* participating in an anomaly 
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now receives a contribution 

oc c/'^^'^BV. (1) 

So, in addition to the diffusion and Ohmic current, there is a new, dissipationless contribution 
proportional to the magnetic field and a density of charges. It is important to note that this 
modification of the constitutive relation occurs regardless of whether the current itself is 
anomalously non-conserved. For example, in the massless QED mentioned above, both the 
axial vector current and the electromagnetic vector current receive this new term in their 
constitutive equations. This result has been obtained previously, by Alekseev, Cheianov, 
and Frohlich, for the cases of 2D field theory, and massless QED j3|. In this work, we will 
confirm their result, using different methods, as well as extending our analysis to a wider 
variety of systems. 

We use two complimentary approaches to elucidate the impact of quantum anomalies 
on the hydrodynamic regime. The first approach uses the gauge/gravity duality, in which 
the hydrodynamic behavior of currents in 4D thermal gauge theories is obtained from the 
dynamics of the dual Yang-Mills fields on black-brane backgrounds in higher-dimensions. In 
this approach the 4D quantum anomaly has a very simple dual description in the higher- 
dimensional theory — it corresponds to the 5D Chern-Simons term in the gauge action. This 
approach applies only to strongly coupled gauge theories with gravity duals, in particular, 
to A/" = 4 super- Yang-Mills theory in the large A^, large t'Hooft coupling regime. In the 
second approach, which is appropriate in the weak-coupling regime, we use a single parti- 
cle spectrum treatment, similar to the "level-crossing" pictures usually used for explaining 
anomalies [J], to examine the behavior of currents in a weakly coupled abelian gauge theory 
with axial anomaly. In both cases we perform our analysis in the background of a constant, 
homogeneous magnetic field and obtain the same additional term in the constitutive relation 
for an anomalous current. We then argue that the form of the new term in the constitutive 
equation that we found is exact, i.e., independent of the strength of interactions. 

The paper is organized as follows. In the first section, we use simple symmetry arguments 
to motivate the inclusion of terms of the form in |3] in the hydrodynamics of anomalous 
theories. In the second section, we explore the R-charge anomaly in the hydrodynamic 
regime of A/" = 4 SYM using a "membrane paradigm" treatment of the dual 5D theory which 
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is very similar to that used in j5[ to examine non-anomalous hydrodynamics in theories with 
gravitational duals. In the third section we present a quasiparticle analysis of the axial 
anomaly in massless QED at finite T and finite (electron) chemical potential. Next we 
present our argument for the universality of this result. This argument is based on the 
equilibrium form of a charge distribution in a system with non- vanishing gauge field. In 
section four, we consider the application of our results to QCD at high temperature and 
strong magnetic field. Then, we conclude. 

II. SYMMETRY CONSIDERATIONS 

The simplest model where the problem of anomalous hydrodynamics appears is QCD with 
two massless quark fiavors. The conserved currents in the theory are the isospin current 
jafi _ g^fj-z^g^ ^]-^g axial isospin current j"'^ = q'J^'y^^q, and the and the baryon current 
Jb ~ Ql^Q- We turn on a background magnetic field B coupled to the baryon current and 
discuss the hydrodynamic behavior of the isospin and the axial isospin currents. 

In the absence of the magnetic field B and at temperatures higher than the chiral phase 
transition, the hydrodynamic equations for the vector and axial charge densities p° and p^" 
are the diffusion equations which can be written as the conservation laws 

p'^ + V-f = 0, p^" + V-j^'' = 0, (2) 

coupled with the constitutive relations 

f = -DVp", j^" = -DVp^" (3) 

We recall that the form of the constitutive equations are dictated by the symmetries (rota- 
tional, C and P) and by the fact that we limit ourselves, in linearized hydrodynamics, to 
terms linear in fields with the lowest number of spatial derivatives. 

However, in the presence of the external magnetic field B, it is possible to write additional 
linear terms 

j'^ = -DVp" + cBp^^ j^'^ = -DVp^" + c'Bp" (4) 

The equations are obviously rotationally invariant. To see that they respect C and P invari- 
ance one recalls that under C —j^, j^^ —>■ j^^, and B —>■ — B, and under P p —>■ p, 
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j -j, ^ -P^ ^ j^ and B ^ B. 

The coefficients c and c' are not fixed by symmetries. In the next two Sections we compute 
these coefficients in some simple theories. 



III. AN APPROACH FROM GAUGE/GRAVITY DUALITY 
A. Introduction 

In this section we follow an analysis developed in Ref. j5|], which employed the 
gauge/gravity duality and the black- hole "membrane paradigm" to demonstrate hydrody- 
namic behavior in a variety of finite-temperature theories with holographic gravitational 
duals. We will work with the particular case of an = 4 super- Yang-Mills theory, believed 
dual to a stack of black D3 branes in AdS^ x 5*^. (The hydrodynamic behavior of this theory 
has been found directly.) For the sake of completeness we will repeat many of the arguments 
used in Ref. This repetition also serves to demonstrate explicitly where our assumptions 
do or do not differ from those of the previous work. 

The AdS^ metric for our black brane configuration is. 



dsl 



(nTR) 



2 r 



u 



i=l 

Note the use of the dimensionless radial coordinate, u, which is related to the usual radial 
coordinate r and horizon location tq through u = rg/r^. The theory possesses SU{A) con- 
served R-charges, whose correlation functions are computable using AdS/CFT. The the 5D 
bulk action of the gauge field dual to the R-current, complete with Chern-Simons term 
is, 

S=-^J d'x v^F;,F^- - j d'x d'^'^e^^^P'^ Ald.A\d,Al. (6) 

In 5D, the Chern-Simons (CS) term is cubic in fields. Thus, two point correlation functions 
can be calculated without this term, but three and higher n-point functions receive anomalous 
contributions from the CS term. The three-point functions can, in principle, be computed 
from the closed-time-path AdS/CFT prescription. In this paper, we compute two-point 
functions in the presence of a background magnetic field — which couples to the two-point 
functions through the CS term in the 5D action. 



In Ref. Pi the general idea, originated from the black-hole membrane paradigm, begins 
with defining conserved currents in terms of field tensors evaluated on a stretched hori- 
zon. The ^escription for how to define these currents is lifted directly from the membrane 
paradigm 6]. Then, one takes the long wavelength limit of the field equations in order to de- 
rive the constitutive equation for the currents. It was shown in Ref. 5| that the 5D Maxwell 
equations imply Pick's law from boundary currents, with a diffusion constant matching that 
found by direct AdS/CFT calculation in the case of A/" = 4 SYM. 

Working in the classical regime of the higher dimensional theory, we use the abelian field 
strength and ignore the /"'"^ terms in the Yang-Mills action. The modified Maxwell equations 
obtained from our action are, 

: [v^i^-'^l + r-y ^'^'-FlF^. = 0. (7) 



In order to simplify our task, let us turn on a constant, homogeneous background magnetic 
field, B, in the 4D theory and discuss charge diffusion on this background. Turning on B in 
4D means imposing a boundary condition that the 5D field strength approach B as m ^ 0. 
One can see that these equations support a constant magnetic field in the three space-like 
dimensions perpendicular to the branes, 

-\e'''F,, = B\ (8) 

We will neglect the backreaction of this field on the background metric. This is justified 
when B <^T^ . Later on, we will work in As = gauge and make use of a partial Fourier 
decomposition of the into plane waves parallel to the horizon with dependence on the 
fifth coordinate left explicit: 



Throughout, we will use the notation 



B = , „ ; ijj 



2{TxTf ' 27rT ' ^ 27rT' 

In addition, we will make the assumptions cl; ~ ~ 5^ and g ^ 1 which in the end are 
seen to be consistent with the modified diffusion equation we obtain. 
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The solution to the field equations describing a constant homogeneous magnetic field, 
together with the metric (0), defines the 5D classical background about which we linearize 
in small perturbations. The linearized Maxwell equations become: 

7^2 Mbc 

dA^99"'9'"F;^] - ^Q^,^g B\F^^ = 0, (9) 

at2 jabc 

dA^99''9''''F;„] + ^^^^^ B\F^^ = 0, (10) 
dAV^9''9''''F^ - j^-^^B^F^, = 0. (11) 



In a theory with anomaly, the definition of current becomes ambiguous. We shall require 
that all currents be gauge invariant under the abelian subgroup of SU{4:) singled out by the 
magnetic field. Thus, we employ the regularization where anomalies appear only in R-charge 
currents for which no external gauge field is turned on. In this regularization the divergence 
of an anomalous current is [3]: 

/V2 _ 1 

In Ref. j^l the radial Maxwell equation, Q is found to function as a conservation equation 
for the membrane paradigm currents. In our case, we define currents for which the radial 
Maxwell equation serves as the Adler-Bell-Jackiw anomaly equation, ()12j). Ultimately, we 
will show that when we take 



Ja — 2 \u^h^ y^'^) 

9SG 

fa = -^^(^f - d^'"'BlF!%^,. (14) 
9sG 

to define the currents for our system, equation (fT^ is properly satisfied. One can already 
see that when B = these currents are conserved by the Maxwell equation. The second 
term on the right hand side of (|TH) is a modification to Pick's law, arising from the presence 
of the Chern-Simons term in the action. 
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B. Analysis 



We now embark on our study of linearized perturbations about our solution. Our general 
strategy will be to demonstrate the same two essential facts in P]. First, as a direct result 
of incoming wave boundary conditions at the horizon 

F^iiush) oc Ftiiush)- 
Second, in the hydrodynamic regime 

Fti{ush) ~ -diAtiush)- (15) 

Our analysis differs from - and not only by our inclusion of new terms in the equations 
of motion and our accommodation of a background B field. We also must employ different 
boundary conditions on At in order to allow for a non-vanishing E in the boundary theory. 
In j^, — > at the boundary leads, through (fT3j) to Fu ~ —diF^ti which supplies the 
diffusion term. In our case, a non-zero At{u = 0) leads, in the same manner, to diffusion 
and Ohmic terms in the constitutive relation, as well as an E ■ B non-conservation term for 
an anomalous current. 

We make use of three Maxwell equations (depending on whether the free index is spatial, 
temporal, or radial) plus two Bianchi Identities: 

dtF,^t - (1 - ^')djF,-, - d'^'^B'^F,) = 0, (16) 



^5^5* - TT^T^^^ifz^-T^^^-^y + d^'^'B'.F^, = 0, (17) 

1 „ 1 



(27rT)2M(l - ' *^ ^ 

.2\ Tpa} ^ o TTia Q Tjia i jabc fyb rpc 



^.[(1 - - (2.T)Mi-^-) "-^- - 1^^^-^- + ^^^^^ - °' ^''^ 

dtF,, - d,F,t - d,Ft, = 0, (19) 

diFtj - d,Fu - dtFij = 0. (20) 

As in [s^ a single wave equation for Fu can be obtained in the near horizon limit. Com- 
bining (|T^ with (fTHj) and (fTBj) respectively we find 

- d,[{l - u'MF.^ + d,Ft^,)] + - d^'^B\dtF^t = 0, (21) 



{2'KTfu{l-v?) ' '''' {2'KTfu 



d^F,% - (1 - u')d,{d,F,% + d,F^^ + d'^'^BptF^^ = 0. (22) 
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Near horizon, these equations simphfy significantly. We proceed under the assumption that 
all three terms of (fTBj) are of the same degree of singularity as we approach the horizon, and 
check the consistency of this assumption after the fact. This allows us to neglect d^F^^, in 
(I22|) and, passing to momentum space, we find 

F^a ^ _ u')\d,F,) + \d^'^B]d,F^^. (23) 
a; a; 

It then follows that for 1 — m ^ uj'^ /q^ the F^t terms in ()2H) can be omitted. The Bianchi 
identity (pUj) indicates that the Fij term can be omitted as well, and we obtain a wave 
equation for Fu: 

d^Fu - (27rT)\(l - u')d,[{l - u')d,F,^] = 0. (24) 
Near horizon, this is solved by 

Fuiu, t) = [a,(l - M)f + A(l - u)-f ]e'"*. (25) 

If we specify incoming wave boundary conditions at the horizon, only the first term can 
be allowed to contribute, and taking time and radial derivatives of the last equation gives 
us 

dtF,^^-{47iT){l-u)d5Fu = 0. (26) 
Bianchi identity (fTIHl then indicates, through (j^!^ that 

^''^-- (4.TKl-.) l-°- 
Since finite energy solutions must decay with time, we have 

Ff 

[AttI j[l — u) 

Note that our results for Fu and F^i are consistent with the assertion that all terms 
of are comparably divergent as we approach the horizon. We have now established 
a relationship between F^i C ji and Fu. To demonstrate a modified Pick's law, we must 
still relate this to a gradient of F^f This is particularly straight forward in the = 
gauge, as F^t = f^sA- Below, we have separeted Ai into two pieces: Ai = A*(x) + Aj('u,t,x) 
where A* gives rise to the constant magnetic field, and Ai{u,t,x.) is the arbitrarily weak 
perturbation about the classical background. In the hydrodynamic regime, and for a weak 
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enough magnetic field, it is possible to find solutions for At{u) and Ai{u), perturbatively in 
u, q, and B, such that 

^^^'^ ' constant, (29) 



while 

Fulsh ~ -^^At\sh. (30) 

Specifically, we take g ^ 1 and u q"^ B"^. 

Finding the perturbative solutions to first order in the small quantities is no simple matter 
for general d"'^'^, but it will not be necessary for our purposes. We will simply show that the 
leading terms of both solutions satisfy the above conditions, determine the value of the 
constant, and demonstrate that corrections to these leading terms are small enough that the 
proceedure is valid. As in [^| we take the streched horizon to be close enough to the actual 
horizon to satisfy 1 — Ugh ^ 1 without being exponentially close: 

-Ljln{l-u,h) -^l. (31) 

Then, all the way down to the stretched horizon we can take 

{l-u)f^l. (32) 

To isolate dependence on radial coordinate, we write the momentum space Maxwell equa- 
tions, using a prime to denote radial differentiation: 

CO Af + (1 - u^)qjA']' + (r^'^B^iqjAl + c^A^ = 0, (33) 

+T^^ff^(liA' - qiJ^) + d"''=B.Mf = 0. (35) 

Then, solving equations (j34|) and (j35|) for uj = q = B = Q with the boundary conditions 
Aj(0) = , At^ = const. , Afj,{u = 1) = const, we find 

Af\u)=uCf^ +A'^'^'\0), (36) 
{l-u^)Af'^' = Cf\ (37) 
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Thus the constant relating At to d^At is one. Now, if we substitute (PHjl into (keeping 
in mind we need only near horizon results to make statements concerning the currents) we 



can relate Cj^^^ to Ct^^\ 



uAf + (1 - u')q^A']' - i(27rT)rf"''^5j[(l - u')Af] = 0. (38) 

In fact, we have a matrix equation, 

Af = -[qjS"' - i(27rT)rf"''"S^][(l - u^)A''], (39) 

that could be used to find A^^^ from Af'^ — though we need not do so here. It is important 
to our argument that this matrix equation be non-singular. While we are not aware of any 
generally applicable reason it should not be, we can always restrict our analysis to some 
subgroup of the R-charge SU{4:) for which this will be true regardless of the relative value of 
q and B. Thus, keeping to our assumption that a) ~ ~ B'^, we have C"'-°'' ~ ^(j'^i^) g^j^^ 



Af)~Af)-ln(l±^). (40) 
■' q 1 — u 



Now we substitute At = Af''' + At^\ Aj = A^p + A^^^ into equation (p^ . again exploiting 
and find 



Af^" = 2q,Af^' - d''''B''Af^'. (41) 



Hence, 

so, 

A^'-CoAfl (43) 

Finally, invoking (jHSl) again, we see 

A^^A^^H'j^)- (44) 

■' q 1 — u 

Thus, we can safely conclude that Fti\sh ~ ~diAt\sh- 
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C. Summary 



In this way, we find ourselves in possesion of solutions satisfying ()29|) and (jHUj) . We 
now see that defining currents as in (jl4|) allows to function as the continuity equation. 
Application of now results in a modified Pick's Law: 

j'^ = -D Vp" + d E" + (r^^&p". (45) 

Here, D = is the diffusion constant, while x = ^g^^ is the susceptibility and a = 
the conductivity P]. As a result, (fT^ becomes a modified diffusion equation: 

dtp" - DV^p" + aV ■ + ( — )d"^"B'' ■ Vp" = ^B'' ■ E^ (46) 

As noted at the begining of this section, we cannot allow currents coupled to a gauge field 
to be non-conserved lest the theory be inconsistent. Hence, no one of the fifteen R-charges 
in the boundary theory will have all of the above terms in its diffusion equation. Charges 
coupling to a gauge field will have an Ohmic term, but no E ■ B term: 

dtp" - DV^p" + aV ■W + ( — )c/"^^B'' ■ Vp" = 0. (47) 

Anomalous charges will have no Ohmic term: 

dtp" - DV^p" + ( — )rf"''"B^ • Vp" = d"^" -B'' ■ E^ (48) 

lon^X I6n^ 

IV. CONSTITUTIVE EQUATIONS AT WEAK COUPLING 

Above, we demonstrated an anomalous modification of constitutive relations in a strongly 
coupled conformal field theory with a gravitational dual. We now turn to theories at weak 
coupling and argue that the same modification should also appear there, and explain its 
physical origin. 

To keep the discussion simple, let's consider the theory a massless fermion. To have a 
hydrodynamic behavior at finite temperature the fermion should interact with itself, but we 
shall assume the interaction to be arbitrarily weak. In accordance to the discussion in the 
previous section, we turn on a background gauge field coupled to the fermion current, 

= (0,0,5x,0). 

12 



Let us emphasize again that A'* is only a background gauge field; we do not include a 
dynamical U(l) gauge field into the theory. The fermionic Lagrangian for our theory is 



= V^T^^'^V', (49) 

with a four-component Dirac spinor and D'^ = — igA^. The fermion Hamiltonian can be 
decomposed into the left- and right-handed parts (Below, a is the vector of Pauh matrices.) 

H^ = iJ (t/SlO- • DV'L - V^rO- • DV'r). (50) 

Energy eigenstates of definite chirality can now be found which satisfy —ia ■ D^^r — E^ip^ 
and icr ■ D'^l = -E'lV'l- The solutions can be separated via 

^R,L(x,y,.)= ( ''^''^'^Me^('=^^-^'=-). (51) 
\0r,l(x) J 

We have two pairs of coupled first ordered differential equations for (p and (p which can be 
written as a second order differential equation for any one of them. For example, 

[V^ - e^B^x^ + E^- kl + eB](t>R = 0. (52) 

where x — x — ky/eB, Thus, we may write solutions for 0r as, 

/ n N 1/4 

'^'^'^''^=[W^) '^^n{V^x)e-^^\ (53) 
Here, Ti.n{x) is the Hermite polynomial defined by 

Jn 

Unix) = {-lfe^'—e^\ (54) 

The energy spectrum is -ER(n) = -^^Jk1 + 2neB, where n is the Landau level's label. There 
is, however, a subtlety at n = 0. The equation from which ipn must now be determined is, 

[k^ - Eji{n)]^R = i(Vs + e5x)0R. (55) 

When n — this equation can only be solved if -E'r(O) = —kz- Thus, the energy eigenstates 
with n = are chiral: the right-handed excitations form a single branch E = —kz, while the 
left-handed energy levels are E = kz- 
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FIG. 1: A non-zero electric charge density gives rise to opposing fluxes of left and right handed 
particles, culminating in Jr-l oc —fjB. 

In the vacuum, all negative energy states are filled. If one turns on a chemical potential 
> for the vector charge, then in addition to the Dirac sea all energy levels with < E < n 
arc populated. For the n = energy levels, this means left-handed fermions with < fc^ < 
and right-handed fermions with — //, < fc^ < are populated. But as all left-handed fermions 
have positive and right-handed ones have negative k^, the net result is that there is a 
nonzero axial current that comes from the n — states. One can also check that the n 
states do not contribute to the axial current, since these energy levels are not chiral: the 
contribution from left- and right-handed sectors cancel each other. The effect survives at 
finite temperature. 

The effect can be quantified most easily by considering a finite volume. If wc place 
the theory in a square box of side length L, with periodic boundary conditions, then the 
momentum runs discrete values ki — ^rii/L. Since the energy eigenvalues do not depend 
on ky, there is a degeneracy of various ky states at each value of E{n). Specifically, we 
requires x = to be inside the box forces < ky < eBL, which means that each value of n 
corresponds to eBLF'/{2'n) possible values of ky. 

In thermal equilibrium the total current is the sum over all energy levels. 



Th — kz 

As mentioned above, states with n 7^ do not contribute to the axial current, since these 
positive energy states exist in equal number with both signs of kz, regardless of chirality. 
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FIG. 2: A chiral density composed of equal numbers of right-handed particles and left-handed 
antiparticles gives rise to a flux of positive charge in the direction of B. 



The function / is the Fermi-Dirac distribution function, 

fR,L{kz,ky;T,fj.) 



(56) 



e/9[|fc.|-M] + 1' 

After taking into account the degeneracy factor related to ky, summing over kz, and regu- 
larizing by subtracting off the contribution of the Dirac sea, one finds 

e 



J5 — JR — JL 



2n' 



(57) 



which coincides with the result obtained for the gravity dual theory in the previous section. 

The same analysis we just used can also be applied in determining the electric current at 
equilibrium, in the presence of a non-zero chemical for the axial charge /is and homogeneous 
background B. Just as in the case of A/" = 4 SYM, we find that the constitutive relation of 
a non-anomalous current is altered in the same manner as that of the anomalous current. 
Specifically, 

e 

Jem = ---^/isB. (58) 
n 

(See Figure 121) This is the exact same result obtained in p. 



V. COUPLING INDEPENDENCE 

The presence of a new type of constitutive terms in the hydrodynamics of quantum field 
theories, on the background of a homogeneous magnetic field, is to be expected on the basis 
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of simple symmetry arguments. These terms look like, 

ja = CabcJb^c ■ 

c 

In the case of a strongly t'Hooft coupled, large N gauge theory dual to an AdS^ super-gravity, 
we see these terms arise naturally from the R-charge anomaly in a membrane paradigm ap- 
proach. In this strongly coupled case, the coefficients of the new terms are simply determined 
from the anomaly coefficient appearing in the divergence equation, and the susceptibility for 
the charge density in the new term. Specifically, 

indicates that, 

'^ac ^abc ■ 

Xc 

(Since B is only present for one charge, C need not depend on the index c. Note also, that 
while the order of indices in unimportant in the tensor structure of the anomaly, it is none 
the less used in Cabc to indicate that we are considering the correction to the current, a, 
that is proportional to the charge density, b. The factor rj is dimensionless, and geometric in 
origin.) In the case of arbitrarily weakly coupled QED at finite temperature, we see the exact 
same contribution to currents arise from a local thermal equilibrium treatment, consistent 
with P] . Here, the currents arise as a direct result of the effect of the chiral anomaly on the 
spectrum of the non- interacting theory. In this section we will formulate an argument as 
to why the coefficients of the new terms should generally be determined as they are in the 
two, disparate examples we have studied above. We'll do this, by combining linear response 
theory with hydrodynamics, and demanding that the hydrodynamics be consistent with the 
zero frequency, zero momentum limit of a thermalized local quantum field theory. 

We use linear response theory to express currents and densities in terms of retarded two 
point correlators, and perturbing source fields j^. Then, we invoke symmetry arguments 
to build hydrodynamic equations for the currents. The hydrodynamic equations will then 
function as equations for the retarded correlators. Solving these equations and setting the 
energy to zero, we find that enforcing the correct behavior of the retarded two point functions 
in the zero momentum limit constrains the coefficients of the constitutive terms to be exactly 
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as they are in our two specific cases above. Doing that we assume the theory to contain 
no massless excitations that could cause the retarded two point functions to be divergent 
at exactly zero energy and small momentum. Specifically, this will mean that no Goldstone 
modes or massless, dynamical gauge fields can be present. The fact that the zero momentum 
limit enforces this constraint indicates that these "kinetic coefficients" for the new terms are 
actually determined by equilibrium physics, as observed in 

We will work with a toy model containing two global charges, one vector and one axial- 
vector, participating in an anomaly with the same background (vector) magnetic field. Be- 
cause there are only two global currents participating in the anomaly in this system, we can 
simplify our notation via Ca = Caw and C\ = Cvav- The same analysis can be easily 
extended to a system with a more general set of currents. More will be said about this at 
the end of the section. 

LRT gives the response of a current to the presence of a source field, up to linear order 
in the source, as 

Here, 

[nr(^ - y)U = o{x' - y'nfM.fM]). 

is the retarded current-current correlator in the absence of the perturbation, and 

is the difference between the current's expectation value in the presence, and absence of the 
source. Here, we will employ a source for the temporal component of our vector current only. 
For consistency with the rest of the paper, we will use the notation, pa = Sj^ in this section, 
thus: 

Pv = ravvA^ ; pA = mAyAl 
Substituting the LRT charge fluctuations into our modified diffusion equations gives, 

{dt - DyV^)py = -OryV ■ Ey - CyBy ■ VpA, (59) 

{dt - DaV^)pa = -CaBv ■ Vpv + VdAYvBy ■ Ey. (60) 
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Since Ey = —'VA^, passing to momentum space allows us to drop a factor of A'^lu, p) from 
each term in both equations, leaving us with equations for the current-current correlators. 

{zu - Dvp')fflvv = ^^v/ + ^Cvp ■ Bv[n°°]Av, (61) 
{tu - DaP^)[II^t^]av = iCaV ■ Bv[n^°]vv - ^r^c^AvvP ■ By. (62) 

Solving for both correlators while keeping (for consistency) only terms up to linear order in 
the magnetic field, we obtain 

~ iCkV ■ Bv[ap^ - ^ciAvv(^^ - Dav")] 

i^R AY - T. T\ 7\ T\ 5 lOdj 

and, 

P"lvv - J-^y (64) 
We are now equipped to ask what zero momentum behavior the hydrodynamic equations 
instill in the retarded correlators, and to demand that this be consistent with the known 
behavior from finite temperature field theory. First, we will take — > and note that in 
this limit, the retarded correlator becomes an analytic continuation of the euclidean time 
(or Matsubara) two point correlator j^. Now, taking p — we note that limp^oin^'jAv is 
well defined. To examine the momentum dependence of our hydrodynamic [11^^] ay it will 
be convenient to separate the momentum into components perpendicular and parallel to the 
magnetic field. Writing B = |Bv| we find. 

Observe, that if we take p — along a contour that keeps the ratio y- constant, the value 
of limp_^o[n5;°]AY will be entirely dependent upon what the value of ^ is. Thus, the limit is 
not well defined, unless the numerator itself is zero. The only way this can happen in the 
context of our hydrodynamic equations, is for the following equality to hold: 

Oa = —f] "AYY- 

0-Y 

Taking the same limits for n^°]vY returns the susceptibility, as it should, and gives no 
information about Ca or Cy. The constant Cy can be found by repeating the same steps 
using only a non-zero A^^^, and examining H^'^jvA- 
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As mentioned earlier, this analysis follows through for more general sets of vector and 
axial-vector currents as well. The currents are fluxes of global charges, of the form 

where Va and are the generators, in flavor space, of the symmetries giving rise to each 
charge. (We will use early Latin indices for vector current generators, and late Latin indices 
for axial current generators, for clarity in what follows.) In general, we can expect [Va, H] 7^ 0, 
[Ax, Ay] 7^ 0, and [Va, A^] 7^ 0. The applied magnetic field couples only to one of the vector 
charges, which we specify here by the index b. In this case, symmetry will demand that new 
constitutive terms have the forms, 

ix = ^ ] CxbcPb^c , ja = ^ ] CazcPz^c- 

b c 

The equations of motion to be solved will now be matrix equations in the space of the global 
charges denoted by {x,a}. Again, one finds that in order for the infrared limit of the two 
point functions to behave properly, 

Cxbc = —f]dxbcXb ^1 

must hold for each possible value of x, and c, while 

Cazc V^azcXz ' 

holds for each a and z, in the presence of Be. This line of argument applies regardless of the 
coupling strength of the theory under consideration. 



VI. QUARK GLUON PLASMA IN MAGNETIC FIELD 

It would be nice to get a look at the kind of effect our modified hydrodynamics could 
have in a real physical system. One case in which transport behavior in a magnetized system 



with anomaly may be of interest is hot QCD in a strong magnetic field 

We will take a look at chiral QCD with three massless flavors, in the presence of a 
homogeneous background U{1)em magnetic field. This theory possesses an axial anomaly 



19 



coupling the electromagnetic (Q), baryon (6), and chiral EM (5) currents. To be explicit, we 
write these currents in terms of the diagonal generators of SU{3)f: 



(66) 
(67) 

(68) 



Where, 



1 



Q = - 

^ 3 



V 



2 
0-10 
-1 J 



^10 0^ 



1 

1 

There are two more currents diagonal in flavor which may be defined as 



v 



J 



d 



with 



\ 




1 
0-1 

These currents have hydrodynamic behavior, and arc anomalous, but their hydrodynamic 
equations do not couple to the other currents, since Tr [A AQ] = 0. Thus, we ignore them. 
The anomaly structure generates non- vanishing three point functions of the types, (jqJqJ5 ) , 
and {fQj^f,). 

In order that we be able to apply the methods of the last section, it is important that 
we take the system to be above the temperature at which the chiral phase transition takes 
place, so that no pion modes will be present to couple to the chiral current. Furthermore, wc 
must treat the U{1)em coupling as arbitrarily small, so the effects of dynamical gauge fields 
will be negligible. Corrections due dynamic electromagnetic fields would be of order 
Under these conditions, the electromagnetic current is diffusive, rather than ohmic. Leptons 
are considered as absent. Thus, all three currents are carried only by quarks, and will have 
identical diffusion constants. 

If we apply the methods of the last section to this model, and again work only to first order 
in the applied magnetic field, we find the following set of coupled hydrodynamic equations. 
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(To keep notation tidy, we will use B' 



TTT— where B is the standard magnetic field.) 



-^B'.V,o 



2XQ 
9X6 



B'-Vj°. 



(69) 
(70) 
(71) 



We take the electric field to be zero here. Thus all three charges are conserved. Passing to 
momentum space, we find a relatively simple eigenvalue equation for the dispersion relations. 



v 



[tu-Dtq^] 

{tu - D^q^) 



z|p B 



-,2XQp.B' -2|p-B' 



[lUJ 



Pq 

VP5/ 



0. 



Clearly, for modes in which the momentum is perpendicular to the magnetic field, the dis- 
persion will be purely diffusive. This will also be the case, when ^ <ti ^, regardless of 
orientation. Modes with momentum parallel to the magnetic field will exhibit more inter- 
esting behavior. Specifically we find the new (normalized) eigenmodes to be, 



1 



1 



Pi = 

P2 = 
P3 = 

with dispersion relations. 



I + XqM Xq/Xb 
1 



Pb + Pq), 



V^S + ixq/Xb) 
1 

\/3 + ixQ/Xb) 



{Pb + Pq + \/l + {XQ/Xb)p5), 
{-Pb - Pq + ^1 + {Xq/Xb)p5) 



toi = —iDp'^, 

2 2 

UJ2 = -iDp - - 
y 



U3 = -tDp' + -^/l + {xQ/xb)p-B'. 



{XQ/Xb)p ■ B', 



(72) 
(73) 
(74) 

(75) 
(76) 

(77) 



Frequency eigenmodes are a mixture of currents defined in equation ()68p. and the eigenfre- 
quencies now have a real component with the same form as the dispersion of Alfven waves. 
(The asymptotic value of XQ/Xb = 2/9 may be substituted in, if exact numbers are desired.) 
An interesting effect of this result, is that an over density of pb + Pq will give rise to a 
dissipationless flow of this charge along the direction of the magnetic fleld. 
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VII. CONCLUSIONS 



We have revisited the famihar "level crossing" picture of the impact of a chiral anomaly on 
the non-interacting limit of a weakly coupled gauge theory. There, we saw that the currents 
participating in the anomaly are altered in a certain way, in the presence of a homogeneous 
magnetic field: j*^ = — ^B^/i'^. We have examined the hydrodynamics of a strongly coupled 
(and strongly t'Hooft coupled) plasma using the Gauge-Gravity duality as a tool. There, we 
saw the R-symmetry currents of large N, Af = 4 SYM receive a new term in their constitutive 
relations, due to their participation in the anomaly: j" = —j^^^d'^^'^B^p^. 

There is a general symmetry argument for the inclusion of such terms in the constitutive 
relations of any anomalous local quantum field theory. In the absence of massless modes 
coupling to the relevant currents, we have a means of determining the coefficients of these 
terms by demanding consistency between hydrodynamics and linear response theory in the 
zero momentum limit. Thus, we conclude that, in the presence of static magnetic fields, 
the hydrodynamic constitutive relations of an anomalous QFT receive a contribution, of the 
form (i'^^'^B^p'^ with a coefficient of — times a geometric factor that can be read off from an 
anomaly equation. This holds regardless of coupling strength, and so long as their are no 
massless dynamical modes coupling to the currents involved. 

If this LRT treatment can be extended to systems with massless modes, it would be pos- 
sible to consider anomalous contributions to magnetohydrodynamics. This is an interesting 
topic for future work, that may significantly broaden the applicability of the considerations 
raised in this paper. 
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